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This  study  reports  predictions  of  the  power  and  fuel  savings  produced  by  thermoelectric  generators  (TEG) 
placed  in  the  exhaust  stream  of  a  sports  utility  vehicle  (SUV)  and  a  stationary,  compressed-natural- 
gas-fueled  engine  generator  set  (CNG).  Results  are  obtained  for  generators  using  either  commercially- 
available  bismuth  telluride  (Bi2Te3)  or  quantum-well  (QW)  thermoelectric  material.  The  simulated  tests 
are  at  constant  speed  in  the  SUV  case  and  at  constant  AC  power  load  in  the  CNG  case.  The  simulations 
make  use  of  the  capabilities  of  ADVISOR  2002,  the  vehicle  modeling  system,  supplemented  with  code 
to  describe  the  thermoelectric  generator  system.  The  increase  in  power  between  the  QW-  and  Bi2Te3- 
based  generators  was  about  three  times  for  the  SUV  and  seven  times  for  the  CNG  generator  under  the 
same  simulation  conditions.  The  relative  fuel  savings  for  the  SUV  averaged  around  -0.2%  using  Bi2Te3 
and  1.25%  using  QW  generators.  For  the  CNG  case  the  fuel  savings  was  around  0.4%  using  Bi2Te3  and 
around  3%  using  QW  generators.  The  negative  fuel  gains  in  the  SUV  were  caused  by  parasitic  losses. 
The  power  to  transport  the  TEG  system  weight  was  the  dominant  parasitic  loss  for  the  SUV  but  was 
absent  in  the  CNG  generator.  The  lack  of  space  constraint  and  the  absence  of  parasitic  loss  from  the 
TEG  system  weight  in  the  CNG  case  allowed  an  increase  in  the  TEG  system  size  to  generate  more  power. 

©  2010  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

New  thermoelectric  materials  show  the  promise  of  reaching 
significantly  higher  values  of  the  thermoelectric  figure  of  merit, 
Z,  and  thus  higher  efficiencies  and  power  densities.  Some  of  the 
most  promising  materials  are  quantum-well  structured  materials 
[1]  and  skutterudites  [2],  Theoretical  analyses  have  shown  that 
quantum-well  (QW)  structured  materials  could  possibly  increase 
Z  by  an  order  of  magnitude  relative  to  the  bulk  material  [3],  Hi-Z 
technologies  has  successfully  made  one  and  two-couple  devices 
[4,5]  and  shown  ZT  values  (where  T  is  the  average  operating  tem¬ 
perature)  of  up  to  4  [6]  compared  to  most  current  materials  which 
are  around  1.  This  has  the  potential  of  roughly  proportionally 
increasing  the  conversion  efficiency  of  thermoelectric  devices. 

In  our  previous  paper  [7],  we  showed  that  it  is  difficult  to  make 
an  automotive  thermoelectric  exhaust  generator  that  generates 
enough  power  to  overcome  the  rolling  resistance  penalty  associ¬ 
ated  with  the  device.  This  is  especially  true  at  low  vehicle  speeds 
when  the  amount  of  power  in  the  exhaust  stream  is  reduced.  The 
purpose  of  this  paper  is  to  examine  the  impacts  that  new  materials 
may  have  on  the  viability  of  vehicle  exhaust  waste  heat  recovery 
with  particular  emphasis  on  how  new  materials  affect  the  behavior 
of  the  entire  system.  Specifically,  we  examine  the  impact  of  QW 
materials  compared  to  Bi2Te3  materials.  QW  materials  were  chosen 
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as  representative  of  the  potential  of  new  materials,  but  skutteru¬ 
dites  also  show  significant  potential.  For  more  information  on 
skutterudite  modules  and  their  optimization,  the  reader  is  referred 
to  the  work  of  El-Genk  and  co-workers  [8-10]. 

There  has  been  a  significant  amount  of  research  done  examin¬ 
ing  waste  heat  recovery.  Hendricks  [11-15]  has  created  a  model 
of  a  generic  generator  and  subsequently  examined  its  optimization 
under  various  conditions  and  with  varying  materials.  In  general, 
these  studies  focus  on  the  thermoelectric  system  including  the 
heat  exchangers  but  do  not  examine  impacts  on  a  vehicle  system. 
Rowe  and  co-workers  have  concentrated  on  the  modules  them¬ 
selves,  doing  fundamental  work  examining  their  optimization  for 
waste  heat  recovery  [16-20],  Additional  work  in  optimization  of 
modules  for  waste  heat  recovery  are  given  by  Chen  et  al.  [21], 
Crane  and  Jackson  [22],  and  Wu  [23]  and  a  discussion  of  the 
impacts  of  novel  materials  on  modules  are  given  in  [24],  The  work 
of  Crane  [25]  includes  some  components  of  a  vehicle  system  by 
accounting  for  power  losses  due  to  the  fan  and  fluid  pump  required 
for  the  coolant  system,  but  this  study  was  done  only  using 
Bismuth-Telluride  modules. 

In  the  following,  we  try  to  assess  the  impact  of  new  materials  on 
the  entire  vehicle  system  including  the  coolant  system,  exhaust 
back  pressure  effect,  and  parasitic  loss  due  to  weight.  Fairbanks 
[26]  gave  an  overview  of  the  application  and  manufacturing  of  QW 
thermoelectric  materials  to  vehicles,  and  Willigan  [27]  discussed 
the  fabrication  of  QW  thermoelectric  material  for  vehicle  waste 
heat  recovery.  In  this  work,  we  investigate  the  potential  benefits 
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Nomenclature 

A 

heat  transfer  area  or  cross-sectional  area  of  thermocou- 

h 

efficiency 

pie  leg 

y 

leg  aspect  ratio 

Ci,C2 

coolant  loss  coefficient  parameters 

p 

rolling  resistance 

Cp 

specific  heat  at  constant  pressure 

kinematic  viscosity 

1 

total  current  of  the  generator 

p 

electrical  resistivity  or  density 

j 

Colburn  factor 

K 

thermal  conductance 

Subscript 

k 

thermal  conductivity  or  pressure  loss  coefficient 

0 

baseline  case 

1 

length  of  thermocouple  leg 

B 

blow-down 

LHVaf 

lower  heating  value  of  the  air-fuel  mixture 

c 

cold  junction 

m 

mass  flow  rate 

D 

driveline 

Nc 

number  of  couples 

E 

engine 

NP 

number  of  unit  generators  in  parallel 

e 

exhaust  side 

Ns 

number  of  unit  generators  in  series 

f 

coolant  side  property 

P 

pressure 

G 

alternator  electrical  generation 

Q 

volumetric  flow  rate 

h 

hot  junction 

R 

electrical  resistance 

j 

module  index 

S 

fuel  savings 

L 

load 

T 

temperature 

oc 

open  circuit 

U 

average  heat  transfer  coefficient 

P 

coolant  pumping 

V 

voltage 

T 

TEG 

v 

vehicle  velocity 

W 

weight 

Acronyms 

W 

power 

CNG 

compressed  natural  gas  generator 

z 

figure  of  merit 

PCHX 

pre-cooling  heat  exchanger 

PCU 

power  conditioning  unit 

Greek 

QW 

quantum  well 

« 

Seebeck  coefficient 

SUV 

sports  utility  vehicle 

A 

change  in  quantity 

TEG 

thermoelectric  generator 

of  QW  thermoelectric  materials  and  try  to  understand  how  much 
of  a  benefit  these  materials  will  have  in  two  practical  applications. 
The  paper  compares  predictions  of  the  power  produced  and  fuel 
saved  when  exhaust  energy  is  recovered  using  a  thermoelectric 
generator  (TEG)  composed  of  one  of  two  materials:  QW  or  com¬ 
mercially  available  Bi2Te3.  The  comparison  is  made  for  a  sports 
utility  vehicle  (SUV)  (Case  1 )  and  a  stationary  compressed  natural 
gas  generator  (CNG)  (Case  2).  To  allow  performance  comparisons, 
simulations  were  performed  for  each  case  with  no  TEG  installed 
(baseline  case),  with  a  TEG  installed  that  used  the  commercially 
available  Bi2Te3  modules,  and  with  a  TEG  installed  that  used  mod¬ 
ules  employing  QW  based  thermoelectric  material.  The  paper  be¬ 
gins  with  a  description  of  the  development  of  a  simulation  code 
followed  by  a  description  and  results  for  the  two  cases. 

2.  Simulation  code  development 

The  simulation  code  consists  of  a  model  of  the  TEG  combined 
with  ADVISOR  2002,  which,  when  originally  obtained,  was  free, 
non-commercial,  and  downloadable.  ADVISOR  is  a  MATLAB®  and 
Simulink®-based  vehicle  simulator  written  at  National  Renewable 
Energy  Laboratory  and  documented  in  [28],  ADVISOR  was  used 
to  model  the  needed  vehicle  inputs,  such  as  fuel  use,  exhaust  flow 
rates  and  temperatures,  alternator  efficiency,  etc.  This  was  coupled 
to  a  model  of  the  TEG  system  we  developed  which  has  been  vali¬ 
dated  [29]  by  comparing  it  to  the  test  results  reported  in  [7]  and 
is  described  herein. 

2.J.  Thermoelectric  generator 

The  TEG  model  is  based  on  an  actual  generator  that  was  built 
using  the  HZ20,  Bi2Te3-based,  thermoelectric  module  sold  by 


Hi-Z  Technology,  Inc  (Hi-Z)  [30],  This  generator  was  tested  in  a 
1999  General  Motors  Sierra  pick-up  truck  [7]  and  is  shown  in 
Fig.  1  without  exterior  insulation,  inlet  and  outlet  transition  pieces, 
and  outer  casing.  As  shown  in  the  figure,  the  tested  generator 
incorporated  sixteen  HZ20  modules  connected  in  series.  Eight  of 
the  modules  were  on  each  side  of  a  central  exhaust  heat  exchanger 
made  of  carbon  steel.  Cooling  was  provided  by  a  pair  of  aluminum 
heat  exchangers,  one  on  each  side.  Thin  aluminum  oxide  wafers 
electrically  isolated  the  thermocouple  circuits  in  the  modules  from 
the  heat  exchangers.  The  contact  surfaces  on  each  side  of  these  wa¬ 
fers  were  covered  with  thermal  grease  and  the  entire  assembly  was 
clamped  together  with  a  200-psi  cold  preload  using  the  shown  pre¬ 
load  subassembly. 

To  model  this  device  we  assumed  thermal  symmetry  with  re¬ 
spect  to  the  orthogonal  vertical  and  horizontal  planes  as  shown 
in  Fig.  1,  and  that  extraneous  heat  losses  could  be  neglected  be¬ 
cause  of  the  insulation  and  casing  mentioned  previously.  The  sym¬ 
metry  planes  divide  the  generator  into  four  zones  extending  over 
its  length.  Each  zone  contains  four  modules  and  may  be  divided 
into  four,  equal-sized  sections  containing  a  single  module,  as 
shown  in  Fig.  2. 

To  model  this  section,  we  neglect  axial  conduction  in  the  metal 
and  assume  that  the  hot  and  cold  junction  temperatures,  Thj  and 
Tcj  over  module  j  are  uniform.  A  steady-state  energy  balance  on 
the  coolant  heat-exchanger  section  considering  convection  in  the 
fluid  and  conduction  heat  transfer  to  the  module  surface  gives 

%%  -  CPfMTfJ+1  -^Vcj  -  Tfj)  =  0,  (1) 

where  Dy  is  the  average  overall  heat  transfer  coefficient  for  the 
coolant  heat  exchanger  for  this  section,  As  is  the  coolant  side  heat 
transfer  area  for  this  section,  thf  is  the  mass  flow  rate  of  the  coolant 
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coolant  symmetry 
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coolant  heat  exchanger 
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thermoelectric  module 


exhaust  heat  exchanger  j 
adiabatic  symmetry  plane 
Fig.  2.  Side  view  of  a  module  section. 


boundary 

K 


through  the  section,  cPf.  is  the  specific  heat  of  the  coolant  evaluated 
at  Tfj,  and  Tf  J  is  the  average  bulk  coolant  temperature  for  this  sec¬ 
tion  defined  by 

Tfj  =  l(Tfj  +  TfJ+%).  (2) 

Similarly,  for  the  exhaust  heat  exchanger  segment 


Wj,  the  power  generated  by  the  module,  is 

Wj  =  IV]  =  I(V0Cj  -  JR,),  (5) 

where  I  is  the  total  current  of  the  generator,  Vj  is  the  module’s 
voltage,  Vocj  is  the  module’s  open  circuit  voltage,  and  R,  is  the  elec¬ 
trical  resistance  of  the  module.  Rj  is  given  by 

Rj  =  Ncyp(Tj),  (6) 

where  Nc  is  the  number  of  (series-connected)  couples  in  the  mod¬ 
ule,  p(Tj)  is  the  resistivity  of  a  couple  evaluated  at  the  average  of 
the  junction  temperatures,  and  y  is  f,  the  aspect  ratio  of  the  legs 
of  a  couple,  assumed  uniform  over  the  generator.  /  and  A  are  the 
leg  length  and  cross-sectional  area  respectively.  The  module’s  open 
circuit  voltage,  Vocj,  is 

V„Cj  =  Xj(Thj-Tcj),  (7) 

where  0Cj  is  the  Seebeck  coefficient  of  the  module,  otj  is  given  by 

oCj=Ncoc(fj),  (8) 

where  oc(Tj-)  is  the  Seebeck  coefficient  of  a  couple  evaluated  at  the 

average  of  the  junction  temperatures.  The  generator  current,  I,  is 


Cp'/ej  ~  cPej  Jej-r -  Thj)  =  0,  (3) 

where  the  variables  have  definitions  similar  to  those  in  Eq.  (1). 
Balancing  the  energy  flow  in  and  out  of  the  module  gives 

UeJAe(Tej  -  ThJ)  -  UfjAf(TCj  -  Tfj)  -  Wj  =  0.  (4) 


14 

Rt  +  Rl  ~r”R,  +  j|» 


(9) 


where  Voc  is  the  total  open  circuit  voltage  of  the  generator,  R  ,  is  the 
electrical  resistance  of  the  thermoelectric  generator,  and  RL  is  the 
load  resistance. 
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A  heat  balance  on  the  hot  surface  of  the  module  as  derived  in 
[31]  is 

UejAe(TeJ  -  Thj)  -  Kj(ThJ  -  TCJ )  -  <XjThjI  +  ±1%  =  0,  (10) 

where  Kj  is  the  thermal  conductance  of  the  module.  Kj  is  given  by 


where  fc(T,j  is  the  thermal  conductivity  of  a  couple  evaluated  at  the 
average  of  the  junction  temperatures. 

Finally,  the  power  of  the  thermoelectric  generator  is 

Wt  =  I2Rl.  (12) 

2.2.  Thermoelectric  module  properties 

The  two  types  of  module  models  used  in  the  calculations  were 
based  on  QW  and  Bi2Te3  materials.  The  Bi2Te3  based  thermoelec¬ 
tric  module  properties  were  obtained  from  the  Hi-Z  HZ20  product 
brochure  [30],  The  HZ20  module  consists  of  71  thermocouples  con¬ 
nected  in  series  electrically.  The  module  properties  were  evaluated 
using  voltage,  power,  and  efficiency  curves  based  on  the  analysis 
developed  in  [29],  However,  the  experimental  power  production 
[7]  was  less  than  the  predicted  power  production  for  the  TEG  using 
standard  HZ20  module  data  [30]  and  experimental  module  surface 
temperatures  [7],  The  Seebeck  coefficient  of  the  HZ20  module  was 
reduced  to  improve  agreement  with  the  experimental  power  pro¬ 
duction.  The  relative  change  in  magnitude  of  the  Seebeck  coeffi¬ 
cient  was  about  8%  at  Th  -  Tc  of  50  °C  to  29%  at  300  °C.  This 
adjustment  resulted  in  a  reduction  of  Z  by  a  factor  of  between 
1.2  and  1.7  over  this  temperature  range.  The  studies  of  the  Bi2Te3 
material-powered  generator,  which  we  name  the  HZ20M,  were 
made  with  this  adjustment. 

The  QW  module  properties  were  based  on  p-type  and  n-type 
Seebeck  coefficient,  resistivity,  and  thermal  conductivity  material 
data,  which  was  furnished  by  Hi-Z  [32],  These  properties  were 
given  as  a  function  of  temperature.  Fig.  3  shows  the  resistivity  of 
the  QW  material  as  function  of  temperature  as  compared  to  Bi2Te3. 
The  Seebeck  coefficient  and  thermal  conductivity  were  relatively 
insensitive  to  temperature.  The  p-type  QW  material  had  ara 
1031  pV/°C  and  k  »  0.5  W/(cm  °C)  as  compared  to  199  pV/°C  and 
0.014  W/(cm  °C)  for  the  p-type  Bi2Te3  material.  For  the  n-type 


QW  material,  the  values  were  a  rj  -1104  pV/°C  and  /i  0. 1 3 
W/(cm°C)  as  compared  to  -175.4  pV/°C  and  0.01  W/(cm°C)  for 
Bi2Te3.  Note  however  that  the  HZ20M  properties  were  not  derived 
from  these  material  properties  but  rather  from  actual  module  per¬ 
formance  tests  as  discussed  above.  Fig.  4  shows  a  comparison  of 
the  Z  factors  of  the  two  modules.  As  the  QW  materials  have  not 
been  produced  as  a  module,  to  reduce  the  number  of  variables 
and  to  provide  a  comparison  to  the  HZ20M,  the  number  of  couples 
and  total  area  for  the  QW  module,  the  HZQW,  were  assumed  to  be 
equal  to  that  of  the  HZ20M.  Fig.  4  shows  that  using  the  QW  mod¬ 
ules  results  in  a  factor  of  increase  in  Z  ranging  from  3  to  6.7. 

2.3.  Fluid  properties 

Both  the  exhaust  and  coolant  fluid  properties  were  evaluated  as 
a  function  of  temperature.  The  density  and  specific  heat  at  constant 
pressure  of  the  exhaust  gas  were  found  from  ideal  gas  mixture  the¬ 
ory.  The  relative  amounts  of  pure  components  in  the  exhaust  mix¬ 
ture  were  determined  using  molar  analysis  for  a  given  fuel  type 
and  the  air-fuel  mass  ratio  obtained  from  the  engine  maps  in  the 
ADVISOR.  The  dynamic  viscosity  of  the  gas  was  found  using  an 
equation  for  gas  mixtures  [33],  For  the  coolant  properties  in  Case 
1,  an  engine  coolant  composed  of  50%  by  volume  solution  of  ethyl¬ 
ene  glycol  in  water,  was  used.  Its  properties  were  evaluated  using 
the  Prestone  Antifreeze  coolant  properties  [34],  In  Case  2,  com¬ 
pressed  liquid  water  at  atmospheric  pressure  was  the  coolant. 

2.4.  Heat  exchangers  and  heat  transfer  coefficients 

The  overall  heat  transfer  coefficients,  Ue  and  Uf,  include  the 
forced  convection  in  the  heat  exchangers,  conduction  through 
the  fins,  and  the  contact  resistance  between  the  heat  exchangers 
and  the  thermoelectric  modules.  The  exhaust  heat  exchanger  was 
an  offset-strip  fin  design  made  of  carbon  steel.  The  geometry  of 
the  fins  is  described  in  [29],  The  design  was  close  to  that  for  which 
correlations  for  the  Colburn  factor  and  friction  factor  were  reported 
in  [35]  to  agree  with  experimental  data  to  within  20%  for  a  range  of 
fin  geometry.  These  correlations  determine  the  Colburn  factor  and 
friction  factor  coefficients  as  a  function  of  the  specific  geometry  of 
the  heat  exchanger  as  well  as  the  Reynolds  number  of  the  flow. 
Thus,  these  parameters  were  evaluated  as  a  function  of  the 
Reynolds  number  of  the  exhaust  flow. 
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The  coolant  heat  exchanger  consisted  of  six,  parallel,  rectangu¬ 
lar  flow  channels  milled  in  a  block  of  aluminium.  The  Reynolds 
number  of  the  flow  based  on  the  hydraulic  diameter  of  these  slots 
ranged  from  about  1300  at  48.3  km/h  to  about  3000  at  1 12.7  km/h 
[29],  Thus  it  ranged  from  the  laminar  to  transitional  region.  How¬ 
ever,  because  of  the  flow  disturbance  imparted  by  upstream  flow 
devices  such  as  the  elbows,  tees,  fittings,  and  manifolds,  we 
assumed  the  flow  to  be  turbulent.  The  Dittus-Boelter  correlation 
[36]  was  therefore  employed  to  calculate  the  heat  transfer 
coefficient. 

2.5.  Solution  strategy 

Eqs.  (1),  (3)-(5),  and  (10)  combined  with  Eq.  (9),  were  simulta¬ 
neously  solved  by  the  multi-dimensional  Newton-Rhapson 
method  to  find  the  four  unknown  temperatures  in  each  section, 

TCj,  Thj,  TeJ,  Tfj,j  =  1,2 . 16,  and  thus  to  calculate  the  performance 

of  the  generator,  principally  its  power  production.  The  analyses 
were  carried  out  at  matched  load,  Rl  =  Rt  for  both  cases.  The  QW 
cases  were  also  run  at  optimum  leg  aspect  ratio,  y.  The  method 
of  determining  optimum  y  will  be  discussed  in  a  subsequent 
section. 

2.6.  Unit  generator 

The  design  shown  in  Fig.  1  was  taken  as  a  “unit  generator”.  The 
number  of  unit  generators  was  varied  depending  on  the  exhaust 
source.  In  the  fixed  engine-generator  study,  Case  2,  for  example, 
the  number  of  unit  generators  was  increased  because  there  was 
more  exhaust  flow  available  and  there  were  no  particular  space 
and  weight  constraints,  as  there  were  in  the  vehicle  of  Case  1.  In 
Case  2,  NP  identical  parallel  lines,  each  parallel  line  consisting  of 
Ns  unit  generators  in  series,  were  fed  by  the  exhaust  and  coolant 
flows.  As  a  result,  each  parallel  line  received  J  of  the  total  flow. 
This  approach  simplified  the  modeling  and  also  kept  the  work 
grounded  in  a  model  that  agreed  sufficiently  well  with  the  exper¬ 
iment  [29], 

3.  Case  1:  sports-utility  vehicle 

In  this  section,  we  examine  the  impact  of  thermoelectric  power 
generation  on  an  SUV  cruising  at  various  speeds.  We  start  with  a 


description  of  the  system  design  followed  by  a  discussion  about 
the  modeling  approach,  determining  of  optimum  coolant  flow  rate 
and  thermoelectric  leg  aspect  ratio,  and  concluding  with  the  dis¬ 
cussion  of  the  results. 

3.1.  System  design 

Fig.  5  illustrates  the  system  schematic  for  Case  1 .  This  scheme 
was  based  on  the  installation  reported  in  [7]  in  which  the  TEG 
was  located  downstream  of  the  catalytic  converters  as  shown  in 
Fig.  5  and  the  intervening  exhaust  piping  was  insulated.  Vazquez 
et  al.  [37]  and  motor  vehicle  industry  contacts  advised  against  locat¬ 
ing  TEGs  between  the  exhaust  manifold  and  the  catalytic  converters 
because  of  possible  slowing  of  heating  of  catalytic  converters, 
resulting  in  an  increased  pollutant  release.  Furthermore,  the  tem¬ 
perature  limits  of  the  module  also  limit  the  possible  positions  for 
the  TEG.  Both  modules  can  be  placed  just  downstream  of  the  cata¬ 
lytic  converter  without  exceeding  their  temperature  limitations. 

Engine  coolant  was  chosen  to  cool  the  TEG  instead  of  air  cooling 
because  the  larger  heat  transfer  coefficient  available  with  liquid 
cooling  allows  smaller  heat  exchangers  to  be  used  and  does  not  re¬ 
quire  fans  at  low  or  zero  vehicle  speed.  Because  the  TEG  is  tapped 
into  the  coolant  system  through  the  cabin  heat  hoses,  the  temper¬ 
ature  of  the  coolant  supplied  to  the  TEG  is  actually  the  highest  in 
the  coolant  system.  The  most  desirable  point  to  tap  into  the  coolant 
system  is  just  after  the  radiator  because  this  is  the  coldest  point  in 
the  system.  However,  the  operation  of  the  by-pass  and  thermostat 
valves  could  then  prevent  coolant  from  flowing  to  the  TEG.  In  the 
experimental  testing  [7],  a  pre-cooling  heat  exchanger  (PCHX) 
was  placed  upstream  of  the  TEG.  Its  purpose  was  to  reduce  the 
TEG  coolant  inlet  temperature  to  levels  more  representative  of 
the  radiator  outlet.  The  effect  of  the  PCHX  on  the  coolant  temper¬ 
ature  drop  was  modeled  by  determining  the  temperature  drop 
across  the  PCHX  at  different  vehicle  speeds  from  the  test  results 
[7],  This  temperature  drop  was  deducted  from  the  engine  coolant 
temperature,  from  ADVISOR,  to  determine  the  coolant  inlet  tem¬ 
perature  to  the  TEG. 

As  shown  in  Fig.  5,  the  electrical  output  of  the  TEG  was  con¬ 
nected  to  the  vehicle’s  electrical  bus  through  a  power  conditioning 
unit  (PCU).  The  PCU  was  essentially  a  direct-current-to-direct-cur- 
rent  converter  functioning  as  a  buck  regulator,  matched  the  gener¬ 
ator’s  output  voltage  to  the  approximately  14.5  V  potential  of  the 
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vehicle’s  direct  current  bus  and  kept  the  TEG  operating  at  its  max¬ 
imum  power  point.  A  PCU  efficiency  of  88%  [38]  was  used  in  the 
simulation.  The  alternator  efficiency  was  50%  [39],  typical  of 
belt-driven  units. 

3.2.  Modeling  approach 

Fig.  6  illustrates  the  data  flow  for  the  system  model.  The  system 
model  consisted  of  three  main  components:  (1 )  ADVISOR  to  model 


the  vehicle  system,  (2)  the  TEG  system  model,  previously  de¬ 
scribed,  and  (3)  the  coolant  system  model.  ADVISOR  did  not  con¬ 
tain  a  model  for  the  engine  coolant  system,  therefore  a  model  of 
it  was  formulated  (discussed  subsequently).  The  basic  input 
parameters  to  the  TEG  model,  exhaust  flow  rate  and  inlet  temper¬ 
ature  and  coolant  flow  rate  and  inlet  temperature,  were  obtained 
from  ADVISOR  and  the  coolant  system  models.  The  output  param¬ 
eters,  electrical  power  and  associated  parasitic  losses,  were  fed 
back  to  ADVISOR  which  then  predicted  the  fuel  consumption. 


thermostat  condition 


Coolant  system  model 


electrical  power 


© 


© 


:k  pressure  from  exhaust  heat  exchanger  & 


transport  TEG  weight 
©additions,  coolant  pump 


ling  approach  Case  1. 
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3.2 A.  Engine  model 

The  engine  model  chosen  from  ADVISOR’S  library  was  a  1991 
Dodge  Caravan  3.0  L  (102  kW)  SI  engine  with  a  torque  scaling  fac¬ 
tor  of  1 .95  applied  to  yield  an  exhaust  temperature  and  mass  flow 
rate  similar  to  those  of  the  test  vehicle  [7],  The  coolant  thermostat 
set  point  for  this  engine  was  modified  from  the  default  value  in 
ADVISOR  to  give  coolant  temperatures  similar  to  those  measured 
during  testing  [7], 

3.2.2.  Coolant  system  model 

This  model  predicted:  the  coolant  flow  rate  through  the  TEG. 
the  additional  pumping  power  needed  to  cool  the  TEG,  and  the 
change  in  flow  rates  through  the  other  components  of  the  coolant 
system.  The  model  was  based  on  the  coolant  system  in  the  truck 
used  for  testing  the  TEG  [7],  The  measurements  taken  during  those 
experiments  and  additional  data  supplied  by  General  Motors  in 
[40,41  ]  were  used  in  the  model  development  and  the  verification 
reported  by  [29],  Fig.  7  is  a  flow  diagram  of  the  system  with  resis¬ 
tor  symbols  representing  loss  coefficients.  The  surge  tank  is  not 
shown.  The  main  flow  loop  is  through  the  pump  and  engine  and 
then  through  four  parallel  paths  consisting  of  the  radiator  by-pass, 
the  radiator,  the  heater,  and  the  TEG.  The  model  is  bi-polar  in  that 
it  treats  only  the  case  when  the  thermostat  valve  is  fully  open  or 
fully  closed.  Sufficient  data  to  model  the  intermediate  condition 
could  not  be  obtained.  However,  the  former  condition  will  exist 
during  cruising  at  constant  speed,  the  state  adopted  in  Case  1. 

The  coolant  was  assumed  incompressible  and  changes  in  poten¬ 
tial  and  kinetic  energy  were  neglected.  The  pressure  losses  for  any 
fluid  element,  Ap,  was  assumed  to  obey 


dictions  of  the  system  model  agreed  well  with  the  available 
measurements. 

GM  provided  data  relating  engine  speed  to  flow  and  pressure 
drop  for  various  components  of  the  unmodified  system,  but  a 
pump  characteristic  curve  at  any  given  engine  speed  was  not  avail¬ 
able.  Therefore,  because  the  flow  rates  through  the  TEG  were  low 
relative  to  the  flow  rates  through  the  other  components  of  the 
coolant  system,  we  assumed  the  variation  in  the  pump  head 
caused  by  the  additional  flow  to  the  TEG  could  be  neglected.  The 
pressure  drop  across  the  TEG  and  PCHX,  A pT  and  the  flow  rate 
through  the  TEG,  Qr,  could  then  be  found  using 

Ap?  =  PfkTQ.T  =  Ap  -  PfliEQ2,  (15) 

where  kE  and  kT  are  the  loss  coefficient  for  the  flow  through  the  en¬ 
gine  and  TEG  loop  respectively.  Ap  is  the  total  pressure  head  pro¬ 
vided  by  the  pump  and  Q  is  the  total  flow  rate.  To  determine  Q, 
the  loss  coefficients  of  the  system  were  used  to  determine  an  equiv¬ 
alent  system  loss  coefficient  using  standard  methods  [42]  and  then 
Eq.  (13)  was  used.  This  process  was  performed  iteratively  to  allow 
for  the  variation  of  the  loss  coefficients  with  Reynolds  number. 

3.3.  Parasitic  power  losses 

The  parasitic  engine  shaft  power  changes  caused  by  the  TEG 
were  due  to  coolant  pumping,  exhaust  blow-down  work,  and  roll¬ 
ing  resistance. 

3.3.1.  Coolant  pumping  power 

The  parasitic  coolant  pumping  power  increase  is  given  by 


A p=pfkQ2,  (13) 

where  Q  is  the  volumetric  flow  rate  through  the  element  and  pf  is 
the  coolant  density.  The  loss  coefficients,  k,  were  assumed  to  de¬ 
pend  only  on  the  Reynolds  number.  Because  it  was  difficult  to  pre¬ 
cisely  define  the  hydraulic  diameter  and  thus  the  Reynolds  number 
in  some  locations,  the  data  was  fit  using  the  following  alternative 
form,  where  v  is  the  kinematic  viscosity, 

k  =  C^fjC\  (14) 

Karri  [29],  reports  the  various  values  of  Ci  and  C2  used  to  fit  the 
experimental  data  for  each  flow  resistance.  Reference  [29]  showed 
that  Eq.  (14)  worked  well  for  most  flow  elements  except  in  the  en¬ 
gine  case,  where  the  scatter  in  the  data  was  large.  However,  the  pre¬ 


AWP  =  WP-  Wpo, 


where  AVVP  is  the  parasitic  pumping  power,  Wpo  is  the  pumping 
power  in  the  baseline  case  (no  TEG),  and  WP  is  the  pumping  power 
with  the  TEG.  For  Case  1,  as  a  full  description  of  the  coolant  system 
was  available,  this  expression  was  used  to  directly  evaluate  the 
additional  coolant  pumping  power.  For  Case  2,  the  coolant  system 
was  generically  defined  and  the  parasitic  pumping  power  was  esti¬ 
mated  by  dividing  the  hydraulic  power  associated  with  the  TEG  cir¬ 
cuit  by  the  pump  efficiency 


pfkr(£ 

W  ' 


(16) 


where  riP  is  the  pump  efficiency.  A  comparison  of  these  two  ap¬ 
proaches  for  Case  1  showed  that  tends  to  predict  a  parasitic 

pumping  power  about  10.5%  greater  than  that  of  (WP  -  VVra). 


Fig.  7.  Coolant  system. 
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3.3.2.  Blow-down  power 

The  blow-down  power  is  the  power  required  for  the  engine  to 
drive  the  exhaust  through  the  exhaust  system.  The  magnitude  of 
this  power  changes  because  of  the  flow  resistance  introduced  by 
the  TEG  exhaust  heat  exchanger.  The  blow-down  power  change 
caused  by  the  altered  exhaust  back-pressure  can  be  calculated  as 


the  wheels  and  is  provided  by  ADVISOR.  The  driveline  transmission 
efficiency  depends  on  vehicle  speed  but  is  roughly  constant  at  about 
0.9.  ADVISOR  also  provides  a  rolling  resistance  coefficient  of  0.012, 
and  the  mass  of  TEG  used  in  the  testing  was  40  kg. 

3.4.  Optimizing  TEG  coolant  flow  rate 


AWb  =  ke  ^f-  keo  — , 


(17) 


where  me  is  the  exhaust  mass  flow  rate,  pe  is  the  exhaust  density, 
and  ke  is  a  lumped  exhaust  loss  coefficient  which  includes  the  losses 
from  the  exhaust  pipes,  muffler,  catalytic  converter,  and  TEG.  A  neg¬ 
ative  or  positive  change  in  blow-down  work  will  depend  on  the  rel¬ 
ative  magnitude  of  the  changes  in  me  and  ke.  If  there  is  a  decrease  in 
shaft  power  due  to  the  gains  produced  by  the  TEG,  me  will  decrease 
and  the  blow-down  work  may  decrease  in  spite  of  the  increased  ke. 
This  was  observed  during  the  testing  reported  by  [7], 

In  the  baseline  case,  the  effect  of  blow-down  power  is  implicitly 
incorporated  through  the  engine  efficiency  map  used  in  ADVISOR. 
To  include  the  additional  effect  of  the  TEG,  the  power  required  to 
force  the  flow  through  the  TEG  is  treated  as  if  it  is  supplied  by  an 
external  fan.  The  power  for  this  fan  is  assumed  to  be  supplied 
directly  from  the  engine  shaft.  In  this  case,  the  additional  blow¬ 
down  work  for  the  TEG  is  given  by 


WBT  =  keT- 1,  (18) 

P2e 

where  keJ  is  the  exhaust  loss  coefficient  determined  from  the  heat 
exchanger  friction  factor  and  we  have  assumed  that  the  (purely 
conceptual)  fan  used  to  supply  this  power  is  100%  efficient. 


3.3.3.  Power  to  transport  TEG  system  weight 

The  parasitic  power  to  transport  the  total  weight  of  the  TEG  sys¬ 
tem  is  given  by 

AVV|)  ^WtV,  (19) 

Bd 

where  pD  is  the  rolling  resistance  coefficient,  Wj-is  the  weight  of  the 
TEG  system,  v  is  the  velocity  of  the  vehicle,  and  rjD  is  the  driveline 
transmission  efficiency.  The  driveline  transmission  efficiency  ac¬ 
counts  for  the  losses  in  the  linkages  between  the  drive  shaft  and 


Increasing  the  coolant  flow  rate,  for  fixed  exhaust  and  coolant 
inlet  temperatures,  exhaust  mass  flow  rate,  and  thermoelectric 
module  properties,  decreases  the  cold  junction  temperature  conse¬ 
quently  increasing  the  junction  temperature  difference  and  the 
power  produced  by  the  modules.  But  an  increase  in  coolant  flow 
rate  will  also  increase  the  pumping  power  proportional  to  the  cube 
of  the  flow  rate,  as  shown  by  Eq.  (16).  The  optimum  flow  rate  is 
that  which  gives  the  largest  difference  between  the  TEG  output 
power  and  the  pumping  power  required  to  circulate  coolant 
through  the  TEG  at  a  particular  engine  rotational  speed.  In  the 
experiment,  the  coolant  flow  rate  was  set  by  opening  or  closing 
butterfly  valves  in  the  TEG  coolant  loop.  In  this  section,  we  theoret¬ 
ically  determine  the  optimal  coolant  flow  rate  (and  thus  the  kT 
value)  for  maximizing  the  net  power  from  the  TEG. 

Studies  were  performed  to  estimate  the  optimum  coolant  flow 
rates  at  various  vehicle  speeds  under  fully  open  and  closed  ther¬ 
mostat  conditions.  Fig.  8  shows  results  obtained  for  the  conditions 
shown  in  Table  1.  The  power  generated  by  the  TEG,  the  power  re¬ 
quired  to  pump  the  coolant  through  the  TEG,  the  net  power  pro¬ 
duced  by  the  TEG,  and  the  change  in  total  flow  rate  and  change 
in  flow  rates  through  the  heater,  by-pass,  and  the  radiator  as  a 
function  of  coolant  flow  through  the  TEG  are  shown.  The  coolant 
flow  rate  that  maximized  the  net  power  produced  by  the  TEG 
was  about  6.3  x  10  5  m3/s.  At  this  flow  rate,  the  changes  in  flow 
rate  through  the  other  components  of  the  coolant  system  were  less 
than  5%.  At  112.7  km/h  the  optimum  was  at  about  1.9  x  10  4  m3/s 
(not  shown).  Based  on  these  studies,  an  intermediate  value  of  kT 
was  chosen  that  resulted  in  a  flow  rate  ranging  from  1.5  x  10  4 
to  2.8  x  10~4  m3/s  over  all  test  speeds. 

3.5.  Optimum  aspect  ratio 

An  increase  in  the  aspect  ratio  of  the  thermoelectric  legs  will 
decrease  the  module  thermal  conductance  and  increase  its  electri¬ 
cal  resistance.  These  competing  effects  produce  a  power  maximum 


TEG  coolant  flow  (m3/s) 


Fig.  8.  Optimal  coolant  flow  analysis  for  Case  1  at  48.3  km/h  vehicle  speed. 
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Data  used  in  the  optimal  coolant  flow  analysis 
shown  in  Fig.  8. 


Variable 


Exhaust  inlet  temperature 
Coolant  inlet  temperature 
Thermostat  condition 


Value 


21.1  °C 
48.3  km/h 
550  W 
0.0142  kg/s 
413.50  °C 

Fully  open 
HZ20M 


Table  2 

Data  used  in  optimizing  yavv  for  Case  1 . 


Vehicle  speed  (km/h)  me  (kg/s)  Te  (°C)  Or  (kg/s)  7)(°C) 

48.3  0.0154  426.4  1.5  x  10  4  79.63 

80.5  0.0286  496.6  2  x  104  77.96 

112.7  0.0520  562.8  2.9  x  1(T4  76.61 


at  a  particular  aspect  ratio.  Consequently,  yow  was  studied  for  both 
cases  to  find  the  ratio  that  maximized  the  total  power  production. 

The  data  used  in  the  optimization  is  given  in  Table  2.  The  ex¬ 
haust  inlet  temperature  and  mass  flow  rate  were  obtained  from 
the  baseline  simulation  results  and  the  coolant  inlet  temperature 
and  flow  rate  were  obtained  from  the  simulation  results  using 
HZ20M  module  properties.  Fig.  9  shows  the  results.  Maximum 
powers  of  117.3,  245.7,  and  459.3  W  corresponding  to  48.3,  80.5, 
and  112.7  km/h  were  obtained  at  the  aspect  ratios  shown.  An  as¬ 
pect  ratio  of  5158  m  1  was  used  for  all  of  the  simulations  using 
the  HZQW  in  Case  1.  The  consequential  reduction  in  the  optimum 
power  at  48.3  and  80.5  km/h  was  about  16.4  and  8.3  W  (14%  and 
3.4%),  respectively. 

For  comparison,  the  HZ20M  legs  have  a  cross-sectional  area,  A, 
of  0.248  cm2  (0.498  cm  x  0.498  cm)  and  are  0.267  cm  long,  giving 
)’hz2om=  107.6  m_1.  Assuming  cross-sectional  area,  A,  for  power- 
optimized  HZQW  legs  equal  to  that  of  HZ20M  would  require  a 
leg  length  of  12.8  cm.  An  alternative  design  procedure  would  be 
to  fix  the  volume  of  the  thermoelectric  material  equal  to  that  of 
the  HZ20M.  Requiring  the  same  leg  volume  as  in  HZ20M  gives 


Aaw  =  .  I^,Ahz2om  =  \/EJo.248  =  0.036  cm2 .  (20) 

y  7q w  »  31-38 

The  legs  of  the  power-optimized  HZQW  module  would  then  be 
1.83  cm  long  with  a  cross-sectional  area,  A,  of  about  0.036  cm2 
(0.189  cm  x  0.189  cm).  (Note  that  Hi-Z,  is  in  the  process  of  building 
QW  modules  having  49  couples,  with  a  leg  length  of  about  1.27  cm 
and  a  cross-sectional  area  of  about  0.258  cm2  [43]  giving  y  = 
492  m  \) 

3.6.  Results 

3.6.1.  TEG  power 

The  power  produced  by  the  TEG  at  vehicle  speeds  of  48.3,  80.6, 
and  112.7  km/h  is  shown  in  Fig.  10.  In  this  figure  and  in  the  next 
two  figures,  the  vertical  ranges  show  the  effect  of  electrical  load 
on  the  system.  At  each  of  the  speeds,  simulations  corresponding 
to  electrical  loads  of  550,  900,  and  1250  W  were  performed.  This 
effect  was  small  relative  to  the  variation  in  power  with  vehicle 
speed.  The  power  increased  with  speed  mainly  because  the  ex¬ 
haust  mass  flow  rate  and  inlet  temperature  increased  with  speed 
while  the  coolant  inlet  temperature  was  relatively  constant.  The 
power  extracted  using  HZQW  modules  was  approximately  three 
times  the  power  extracted  using  HZ20M  modules. 

3.6.2.  Fuel  savings 

An  theoretical  expression  for  the  relative  fuel  savings  will  be 
useful  in  discussing  the  results.  Assuming  that  the  engine  exhaust 
contains  no  liquid  water,  an  energy  balance  on  the  engine  gives: 

WE  =  t]Emc.LHVAF.  (21) 

where  WE  is  the  engine  shaft  power,  rjn  is  the  engine  shaft  effi¬ 
ciency,  and  LHVap  is  the  lower  heating  value  of  the  air-fuel  mixture 
at  inlet  conditions.  The  engine  shaft  power  is  used  (in  general)  for 
driveline  propulsion,  coolant  pumping,  blow-down  power,  alterna¬ 
tor  electric  power  generation  WG,  as  well  as  other  possible  loads 
from  the  air-conditioning  unit  and  the  cooling  fan  which  we  assume 
constant  for  this  analysis. 

WE  =  WD  +  Wp  +  WB  +  Wc+  Wother.  (22) 

The  portion  of  the  shaft  power  required  to  drive  the  alternator  is 


YOtT1) 


Fig.  9.  SUV  yaw  optimizatic 
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Vehicle  speed  (km/h) 

Fig.  10.  TEG  power  production  Case  1. 


WG 


WL 


(23) 


where  WL  is  the  electrical  load  and  rjc  is  the  efficiency  of  the  alter¬ 
nator.  Using  Eqs.  (21)-(23)  and  taking  the  difference  between  the 
baseline  case  and  the  case  with  the  TEG  gives  an  equation  for  the 
fuel  savings  ratio,  S, 


Ai he  _  AWbA-'&W^  ripgjWj 

rhe  ~  w,n  +  nGWm  ’ 


where  rjpcu  is  the  PCU  efficiency  and  we  have  assumed  that  the  elec¬ 
trical  load  remains  constant  but  is  offset  by  the  power  supplied  by 
the  TEG.  Eq.  (24)  shows  that  the  fuel  savings  depend  on  three 
dimensionless  ratios:  (1 )  the  loss  ratio,  which  is  the  ratio  of  the  total 
parasitic  loss  to  the  baseline  shaft  power,  (2)  the  efficiency  ratio, 
which  is  the  ratio  of  the  PCU  efficiency  to  the  alternator  efficiency, 
and  (3)  the  power  ratio,  which  is  the  ratio  of  the  TEG  power  to  the 
baseline  shaft  power.  We  will  use  this  equation  to  help  explain  the 
fuel  savings. 


Fig.  1 1  shows  fuel  savings.  Positive  fuel  savings  were  attained  at 
all  the  vehicle  speeds  using  HZQW  module.  The  maximum  fuel 
savings  of  about  1.5%  was  obtained  at  a  vehicle  speed  of 
112.7  km/h  using  this  module.  However,  negative  fuel  gains  were 
registered  at  vehicle  speeds  of  48.3  and  80.6  km/h  using  the 
HZ20M  module.  Using  Eq.  (24)  and  assuming  no  parasitic  losses, 
for  the  HZQW -TEG  power  varying  from  100  to  450  W  the  relative 
fuel  savings  should  be  about  2-2.3%  and  for  the  HZ20M-TEG  power 
varying  from  33  to  140  W  the  relative  fuel  savings  should  be  about 
0.68  to  0.83%.  However,  as  shown  in  Fig.  1 1  under  both  cases  the 
relative  fuel  savings  is  significantly  lower  because  of  parasitic  ef¬ 
fects  of  the  TEG  on  the  vehicle  system. 

A  comparison  of  the  TEG  shaft  power  to  its  associated  parasitic 
losses  powers  is  shown  in  Fig.  12.  To  base  the  comparison  to  a  stan¬ 
dard  point  of  reference,  the  TEG  power  was  converted  into  its 
respective  shaft  powers  as  shown  in  Eq.  (24).  The  parasitic  blow¬ 
down  power  does  not  include  the  entire  exhaust  system  because 
loss  coefficients,  for  example  for  the  exhaust  pipe  and  muffler, 
were  not  available.  Only  the  effect  of  the  flow  resistance  from 
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Fig.  12.  Parasitic  losses  Case  1. 


the  TEG  is  included  as  shown  in  Eq.  (18).  Fig.  12  shows  that  the  to¬ 
tal  parasitic  loss  exceeded  the  power  produced  by  the  TEG  at  lower 
vehicle  speeds  for  the  HZ20M.  This  resulted  in  a  fuel  consumption 
increase  compared  to  the  base  case,  as  shown  in  Fig.  1 1 .  The  dom¬ 
inant  parasitic  loss  is  the  power  to  transport  the  TEG  system 
weight  which  increases  with  vehicle  speed.  The  difference  in 
weight  for  the  unit  TEG  due  to  the  change  in  modules  from 
HZ20M  to  HZQW  will  depend  on  the  leg  geometry  of  the  QW  mod¬ 
ule.  For  example,  the  power-optimized  QW  module  with  a  leg 
cross-sectional  area,  A(lw,  of  about  0.036  cm2  would  have  a  weight 
of  about  24  g  [44],  We  do  not  know  the  feasibility  of  these  designs 
from  a  manufacturing  standpoint.  In  contrast,  the  weight  of  the 
HZ20M  module  is  about  115  g  [30]  and  the  total  weight  of  the 
16  HZ20  modules,  used  in  the  tested  TEG,  was  about  only  4.6%  of 
the  total  TEG  weight  which  was  dominated  by  the  heat  exchangers. 
This  implies  that  the  change  in  thermoelectric  material  will  not  re¬ 
sult  in  a  significant  change  in  the  parasitic  loss  associated  with  the 
TEG  weight. 

The  coolant  pumping  power  and  the  blow-down  power  are  less 
significant.  The  coolant  pumping  power  is  almost  independent  of 
the  vehicle  speed  because  the  coolant  system  by-pass  valve  main¬ 
tains  an  approximately  constant  pressure  drop  across  the  heater 
core  even  as  the  engine  rpm  changes  [41],  Therefore  the  resulting 
flow  and  pressure  drop  through  the  TEG  coolant  loop  are  approxi¬ 
mately  independent  of  vehicle  speed.  The  coolant  pumping  power 
was  the  same  for  the  two  modules.  The  blow-down  power  is  nearly 
the  same  between  the  two  cases  as  well,  primarily  because  there 
was  not  a  large  change  in  exhaust  mass  flow  rate.  For  the  cases 
in  which  the  TEG  produces  a  net  positive  power,  what  is  shown 
is  actually  a  worst-case  representation.  The  blow-down  power  in 
the  remainder  of  the  exhaust  system  decreases  relative  to  the  base¬ 
line  because  of  the  reduced  exhaust  gas  flow  rates. 

3.6.3.  Energy  budget 

Table  3  compares  the  energy  budget  for  the  HZQW  and  HZ20M. 
The  following  are  evaluated  in  the  table: 

•  Maximum  possible  exhaust  power  that  can  be  extracted, 

We,mra  =  me(cPe  Te.  -  cPefj Tu),  (25) 

where  Wemax  is  the  total  possible  exhaust  power  that  can  be  ex¬ 
tracted,  cPel.  is  the  specific  heat  of  the  exhaust  evaluated  at  Tc-i,  Tei 


Table  3 

Energy  budget  Case  1. 


Vehicle  speed  (km/h)  48.3 


HZQW 

We,max 

We 

Effectiveness 
TEG  power 

Heat  rejected  to  coolant 

TEG  efficiency 

HZ20M 

We,max 

we 

Effectiveness 
TEG  power 

Heat  rejected  to  coolant 
TEG  efficiency 


(kW)  6.27 

(kW)  2.26 

{%)  36.12 

(kW)  0.100 

(kW)  2.16 

(%)  4.41 


(kW)  6.37 

(kW)  2.59 

(%)  40.80 

(kW)  0.033 

(kW)  2.56 

(%)  1.26 


80.5 


112.7 


30.46 

5.31 

17.43 

0.453 

4.86 

8.53 

31.01 

6.21 

20.04 

0.140 

6.07 

2.25 


is  the  exhaust  inlet  temperature  into  the  TEG,  cPe  is  the  specific 
heat  of  the  exhaust  evaluated  at  T/j,  and  Tfi  is  the  coolant  inlet  tem¬ 
perature  into  the  TEG. 

•  Actual  exhaust  energy  extracted, 

We  =  ihe(cP'tTe, t  -  cPl  Je,„),  (26) 

where  We  is  the  exhaust  energy  extracted,  cPeo  is  the  specific  heat  of 
the  exhaust  evaluated  at  Te,0,  and  Te,0  is  the  exhaust  outlet  temper¬ 
ature  from  the  TEG. 

•  The  effectiveness  of  TEG  is  the  ratio  of  exhaust  energy  extracted 
to  maximum  possible  exhaust  energy  that  can  be  extracted. 

•  Power  generated  by  the  TEG. 

•  Exhaust  energy  rejected  to  coolant. 

•  TEG  efficiency,  the  ratio  of  power  generated  by  the  TEG  to 
exhaust  energy  extracted. 

Table  3  shows  the  TEG  efficiency  increased  by  a  factor  of  about 
3.2-3.8  because  of  the  change  from  the  HZ20M  to  the  HZQW.  Thus, 
because  the  HZQW  rejects  less  heat,  a  lower  coolant  flow  rate 
could  be  used.  This  would  reduce  the  parasitic  pumping  losses 
and  further  increase  the  gains  associated  with  switching  to  QW 
material. 

Table  3  shows  the  effectiveness  of  the  unit  TEG  using  HZ20M 
modules  is  about  41%,  29%,  and  20%  at  48.3,  80.5,  and  112.7  km/h 
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respectively.  This  indicates  that  the  exhaust  energy  at  higher 
speeds  was  underutilized.  The  effectiveness  of  the  HZQW  de¬ 
creases  to  36%,  25%,  and  17%  at  48.3,  80.5,  and  112.7  km/h  respec¬ 
tively.  This  highlights  the  fact  that  the  switch  to  better 
thermoelectric  materials  needs  to  be  augmented  with  better  heat 
exchanger  designs.  A  possible  scenario  would  be  to  increase  the 
plan  area  of  the  TEG  system,  so  that  more  modules  could  be 
employed,  thereby  increasing  the  power  generated  and  exhaust 
energy  extracted.  However,  the  size  of  the  TEG  system  was  con¬ 
strained  by  the  limited  space  in  the  underbody.  Additionally,  an  in¬ 
crease  in  system  size  would  increase  the  system  weight  and 
consequently  the  rolling  resistance  parasitic  loss. 

To  summarize,  the  application  of  QW  and  Bi2Te3  based  TEG  to 
the  SUV  has  been  studied.  The  power  generated  by  the  TEG  using 
HZQW  was  about  100-400  W.  This  should  result  in  a  fuel  savings 
of  about  2-2.3%.  However,  the  fuel  savings  decreased  by  about 
30%,  primarily  due  to  the  rolling  resistance  parasitic  loss.  The  ex¬ 
haust  energy  was  underutilized  because  of  space  limitations,  espe¬ 
cially  at  higher  vehicle  speeds. 

4.  Case  2:  natural  gas  engine  generator 

In  this  section,  we  examine  the  impact  of  the  thermoelectric 
power  generation  on  a  stationary,  compressed-natural-gas- 
powered  generator  at  various  electrical  loads. 

4.1.  System  design 

Fig.  13  illustrates  the  system  schematic  for  Case  2.  The  absence 
of  parasitic  loss  from  TEG  system  weight,  the  lack  of  space  con¬ 
straint,  and  the  higher  exhaust  flow,  allowed  us  to  use  more  than 
one  unit  TEG.  The  TEG  stack  was  located  downstream  of  the  cata¬ 
lytic  converter  for  the  reasons  given  earlier. 

As  shown  in  Fig.  13,  the  electrical  output  of  the  TEG  stack  is  con¬ 
nected  to  the  AC  generator’s  electrical  bus  through  a  PCU.  The  PCU 


for  this  study,  besides  matching  the  TEG  stack  output  to  the  AC 
generator  voltage  and  tracking  the  maximum  power  point  of  the 
TEG  stack,  converts  the  TEG  stack  output  from  direct  to  alternating 
current.  The  PCU’s  efficiency  was  assumed  to  be  88%.  The  AC  gen¬ 
erator  performance  data  was  obtained  from  the  specification  sheet 
for  the  John  Deere  Power  Tech  6068H  Diesel  engine  for  generator 
set  applications.  This  AC  generator  was  chosen  based  on  the  load 
capabilities  of  the  CNG  engine.  The  AC  generator  has  an  electrical 
output  of  124-143  kW  operating  at  a  frequency  of  60  Hz  and  an 
efficiency  of  about  88-92%  [45], 

4.2.  Modeling  approach 

Fig.  14  illustrates  the  data  flow  for  the  system  model.  This  mod¬ 
el  consists  of  three  main  components:  (1 )  engine  torque  calculator, 
(2)  CNG  model,  and  (3)  TEG  stack  model.  The  engine  torque  calcu¬ 
lator  determines  the  torque  required  at  the  engine  shaft  for  a  given 
AC  generator  load,  efficiency,  and  shaft  speed.  The  CNG  model  cal¬ 
culates  the  fuel  usage,  exhaust  mass  flow  rate  and  temperature  and 
other  variables  based  on  the  required  torque  and  engine  speed,  and 
the  TEG  stack  model  is  similar  to  the  TEG  system  model  which  was 
previously  described.  Since  ADVISOR  is  designed  as  a  vehicle  anal¬ 
ysis  tool,  changes  were  made  to  the  CNG  model  to  account  for  the 
stationary  characteristics  of  Case  2.  The  engine  used  in  this  case 
study  was  a  John  Deere  8.1  L  (186  kW)  CNG  SI  engine,  the  only 
CNG  engine  available  in  the  library  of  ADVISOR.  The  exhaust  mass 
flow  rate  from  this  engine  was  about  five  times  that  of  the  SUV  in 
Case  1  at  80.5  km/h. 

4.3.  Optimization 

The  optimum  number  and  arrangement  of  TEGs  was  found  by 
calculating  the  net  maximum  power  from  the  TEG  stack  under  var¬ 
ious  combinations  of  parallel  and  series  configurations  using  the 
HZQW.  The  power  from  the  TEG  stack  was  determined  using  an 


,  Case  2. 
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Fig.  14.  Modeling  approach  Case  2. 


optimum  aspect  ratio,  y,  found  for  each  stack  configuration.  Case  2 
losses  consisted  of  coolant  pumping  power  and  blow-down  power. 
The  unit  TEGs  in  series  were  arranged  in  a  counter-flow  arrange¬ 
ment.  The  flow  pattern  for  a  2x2  configuration  is  shown  in 
Fig.  13.  The  exhaust  flow  and  inlet  temperature  conditions  were 
obtained  from  the  baseline  simulations.  The  total  exhaust  mass 
flow  rate  was  0.163  kg/s  and  the  exhaust  inlet  temperature  was 
469  °C.  We  assumed  that  there  was  a  continuous  supply  of  water 
from  city  mains  or  cooling  towers  to  cool  the  TEG  stack.  The  cool¬ 
ant  flow  through  each  of  the  series  lines  and  the  inlet  temperature 
to  the  stack  were  set  at  2  x  10  4  m3/s  and  15  °C  respectively. 

From  the  optimization,  we  selected  an  optimum  stack  configu¬ 
ration  of  five  parallel  lines  with  ten  unit  TEGs  each  in  series.  This 
resulted  in  a  total  stack  power  of  5523,  a  1 10.5  W  of  power  per  unit 
generator.  However,  the  maximum  net  power  produced  by  the 
stack  was  6495.7  W  using  a  total  of  260  unit  TEGs,  at  25  W  of 
power  per  unit  generator.  Although  there  was  no  parasitic  power 
penalty  because  of  this  extra  weight,  it  seemed  unreasonable  to 
employ  TEGs  at  such  a  low  unit  power.  Inclusion  of  factors  such 
as  cost  and  space  constraints  would  have  imposed  further  limits 
on  the  TEG  stack.  However,  the  above  optimization  allowed  us  to 
select  what  seemed  a  reasonable  number  and  arrangement  of  unit 


Table  4 

A  comparison  of  flow  conditions  between  case  studies  1  and  2. 


Variable  Case  2  Case  1 


exhaust  flow  (kg/s)  0.03  0.03 

Exhaust  inlet  temperature  (°C)  468.5  495 

Exhaust  outlet  temperature  (°C)  459  397 

Coolant  flow  (m3/s)  2  x  10  4  2  x  10  4 

Coolant  inlet  temperature  (°C)  28.7  78 

Coolant  outlet  temperature  (°C)  29  81 


TEGs  in  the  stack  for  the  purposes  of  this  analysis.  With  this  num¬ 
ber  of  series  and  parallel  generators,  the  exhaust  and  coolant  flow 
conditions  for  a  unit  TEG  in  a  series  line  are  comparable  to  that  of 
Case  1  at  80.5  km/h  as  shown  in  Table  4. 

The  aspect  ratio  optimization  for  the  TEG  stack  using  five  paral¬ 
lel  lines  with  each  line  consisting  of  10  unit  TEGs  in  series  gave  an 
optimum  of  22,222  nr1.  As  in  Case  1,  Section  3.5,  the  number  of 
couples  and  the  total  area  for  the  QW  module,  HZQW,  were  as¬ 
sumed  to  be  equal  to  those  of  the  HZ20M. 

4.4.  Results 

4.4.1.  TEG  power 

The  power  produced  by  the  TEG  stack  at  AC  generator  loads  of 
124, 130, 137,  and  143  kW  is  shown  in  Fig.  15.  The  error  bars  show 
the  effect  of  AC  generator  efficiency  on  the  system.  At  each  of  the 
AC  generator  loads,  simulations  corresponding  to  AC  generator 
efficiency  of  88%  and  92%  were  performed.  The  TEG  stack  power  in¬ 
creased  slightly  with  AC  generator  load,  because  the  exhaust  mass 
flow  rate  increased  with  AC  generator  load  while  the  coolant  inlet 
temperature  and  mass  flow  rate  were  constant.  The  power  ex¬ 
tracted  using  HZQW  modules  was  approximately  6.5-7.4  times 
the  power  extracted  using  HZ20M  modules. 

4.4.2.  Fuel  savings 

Fig.  16  shows  fuel  savings  for  Case  2.  In  this  figure  and  in  the 
following  figure,  the  vertical  ranges  show  the  effect  of  AC  generator 
efficiency  on  the  system.  The  maximum  fuel  savings  of  about  3% 
was  obtained  with  the  HZQW  module.  This  occurred  at  a  AC  gener¬ 
ator  load  of  124  kW  and  efficiency  of  88%. 

Following  the  generalized  fuel  savings  analysis  in  Case  1,  the 
TEG  stack  power  was  converted  to  its  corresponding  fuel  savings 
using  the  AC  generator  efficiency  of  88%  or  92%  and  a  PCU  efficiency 
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Fig.  15.  TEG  power  production  Case  2. 
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Fig.  16.  Fuel  savings  Case  2. 


of  88%.  The  relative  fuel  savings  for  the  TEG  power  (HZQW)  varying 
from  5.3  to  6  kW  is  about  3.2-3.3%  and  for  the  TEG  power  (HZ20M) 
varying  from  0.74  to  0.9  kW  is  about  0.5%.  These  are  in  good 
agreement  with  the  calculations  based  on  the  fuel  usage  from 
the  simulations  with  and  without  the  TEG,  as  shown  in  Fig.  16.  This 
indicates  that  the  parasitic  losses  were  not  significant,  as  will  be 
verified  subsequently.  Also,  from  Fig.  16  the  fuel  savings,  using 
HZQW  modules,  decreases  with  increasing  AC  generator  load 
because  of  the  decrease  in  power  ratio  as  the  AC  generator  load 
increased.  The  fuel  savings  using  the  HZ20M,  also  decreased  as 
the  AC  generator  load  increased.  However,  this  decrease  was  too 
small  to  appear  in  the  figure. 

The  parasitic  losses  were  evaluated  as  previously  discussed, 
with  that  caused  by  weight  set  to  zero.  The  parasitic  blow-down 
power  is  not  included  because  as  in  Case  1  necessary  data  was 
not  available.  The  effect  of  back-pressure  on  the  shaft  power  was 
accounted  for  as  in  Case  1. 

Fig.  17  shows  a  comparison  of  the  TEG  power  to  the  parasitic 
coolant  pumping  power  and  the  back-pressure  power.  All  the  pow¬ 
ers  have  been  converted  into  their  respective  shaft  powers  as  in 


Case  1.  However,  the  pump  efficiency  was  not  accounted  for  in 
the  evaluation  of  parasitic  pumping  power.  The  pumping  power, 
without  accounting  for  the  pump  efficiency,  was  about  90  W.  A 
pump  efficiency  of  25%  would  result  in  a  parasitic  pumping  power 
of  360  W.  These  powers  correspond  to  a  range  of  10-50%  of  the 
TEG  stack  power  for  the  HZ20M  or  1.6-7%  for  the  HZQW.  The 
blow-down  power  is  again  less,  being  at  the  most  only  8.6%  of 
the  power  generated  by  the  TEG  for  HZ20M  case  or  1.3%  for  the 
HZQW. 

4.4.3.  Energy  budget 

Table  5  compares  the  energy  budget  for  the  HZQW  and  HZ20M. 
The  tabulations  are  the  same  as  those  defined  in  the  energy  budget 
for  Case  1.  From  Table  5,  the  TEG  stack  efficiency  due  to  the  change 
in  modules  from  HZ20M  to  HZQW  increased  by  a  factor  of  about 
10.  A  comparison  of  the  effectiveness  for  the  stack  shows  that,  al¬ 
ready  70-66%  of  the  maximum  possible  energy  has  been  extracted 
using  HZQW  and  this  further  increases  to  96-97%  using  HZ20M 
modules.  However,  a  significant  fraction  of  this  extracted  energy 
is  being  rejected  to  the  coolant,  about  90%  using  HZQW  and  about 
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Table  5 

Energy  budget  Case  2. 


AC  generator  load 
HZQW 


Effectiveness 
TEG  stack  power 
Heat  rejected  to  coolant 
TEG  stack  efficiency 
HZ20M 


Effectiveness 
TEG  stack  power 
Heat  rejected  to  coolant 
TEG  stack  efficiency 


(kW)  124  130  137 


(kW)  80.3 

(kW)  56.3 

M  70.1 

(kW)  5.27 

(kW)  51.0 

(%)  9.36 


(kW)  82.5 

(kW)  80.3 

(%)  97.3 

(kW)  0.72 

(kW)  79.6 

(%)  0.90 


83.5  87.2 

57.5  58.8 

68.8  67.4 

5.46  5.66 

52.0  53.1 

9.50  9.63 


85.7  89.4 

83.2  86.4 

97.0  96.7 

0.77  0.83 

82.4  85.6 

0.93  0.96 


143 


90.2 
59.8 

5.83 

54.0 

9.74 

92.5 

89.2 
96.4 


99%  using  HZ20M  modules.  Any  further  increase  in  system  size, 
though  increasing  the  power  from  the  stack,  might  not  be  feasible 
economically.  In  fact,  with  the  aid  of  better  module  and  heat  ex¬ 
changer  designs  the  system  size  could  be  reduced  while  retaining 
the  stack  efficiency  and  effectiveness. 

To  summarize,  the  TEG  stack  generated  about  53-5.8  kW  using 
optimized  QW  modules.  This  power  was  greater  than  that  of  the 
power  from  Bi2Te3  modules  by  a  factor  of  6.6-73.  The  power  gen¬ 
erated  by  the  QW  modules  resulted  in  a  fuel  savings  of  about  3%. 
This  stationary  application  allowed  a  larger  system  size  and  also 
resulted  in  a  low  total  parasitic  loss.  The  greater  system  size  also 
increased  the  power  generated,  effectiveness,  and  efficiency  of 
the  stack. 

Though  the  change  in  power  generated  between  Cases  1  and  2 
is  quite  large,  the  change  in  relative  fuel  savings  is  not  in  propor¬ 
tion.  As  discussed  earlier  the  magnitude  of  power  ratio,  efficiency 
ratio,  and  the  loss  ratio  determine  the  relative  fuel  savings  for  each 
of  these  cases.  For  Case  1,  the  power  ratio  for  the  unit  TEG  using 
HZQW  modules  at  112.7  km/h  was  about  1.3%,  and  the  efficiency 
ratio,  ticpuMg,  was  1.8.  For  Case  2  with  the  HZQW  generator  at 
143  kW  AC  generator  load,  the  power  ratio  was  greater  at  3.3% 
but  the  efficiency  ratio  was  only  1.0  because  of  the  higher-effi¬ 
ciency  generator  used.  Thus,  although  Case  2  had  a  greater  power 
ratio  and  a  smaller  loss  ratio,  these  impacts  were  mitigated  by  the 


fact  that  the  stationary  electrical  generator  was  more  efficient  than 
the  vehicle  alternator. 


5.  Conclusions 

The  potential  benefits  of  the  application  of  QW  and  Bi2  Te3 
based  thermoelectric  power  generation  from  the  SUV  and  CNG 
engine  power  generator  were  examined.  Under  both  these  applica¬ 
tions  the  QW  based  TEG  generated  more  power  relative  to  the  Bi2 
Te3  based  TEG.  For  the  SUV,  the  QW  based  TEG  generated  about 
100-450  W,  which  should  result  in  a  fuel  savings  of  about 
2-2.3%.  However,  the  fuel  savings  decreased  by  about  30%,  primar¬ 
ily  due  to  the  rolling  resistance  parasitic  loss.  In  contrast,  for  the 
CNG  case  the  total  parasitic  loss  was  low  relative  to  the  TEG  stack 
power.  This  was  primarily  because  this  was  a  stationary  applica¬ 
tion  and  the  lack  of  space  constraint  has  allowed  an  increase  in 
system  size.  The  optimized  QW  based  TEG  stack  generated  about 
53-5.8  kW  resulting  in  a  fuel  savings  of  about  3%.  The  low  relative 
fuel  gains  despite  the  higher  TEG  stack  power  and  lower  parasitic 
loss  was  primarily  due  to  the  high  efficiency  of  the  electrical  gen¬ 
erator.  The  energy  budget  showed  that,  under  both  the  applica¬ 
tions,  a  significant  fraction  of  the  exhaust  energy  was  rejected  to 
the  coolant. 
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